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Abstract. The concept of a generalized nonanalytic expansion which involves nonanalytic 
combinations of exponentials, logarithms and powers of a coupling is introduced and its use 
illustrated in various areas of physics. Dispersion relations for the resonance energies of odd 
anharmonic oscillators are discussed, and higher-order formulas are presented for cubic and 
quartic potentials. 
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1 Introduction and motivation 

In many cases, a simple power series, which may be the result of a Taylor expansion, is not enough 
in order to describe a physical phenomenon. Furthermore, even if a power series expansion 
(e.g., of an energy level in terms of some coupling parameter) is possible, then it may not be 
convergent [1, 2, 3, 4, 5]. Physics is more complicated, and generalizations of the concept of 
a simple Taylor series are called for. 

Let us start with a simple example, an electron bound to a nucleus. It is described to a good 
accuracy by the Dirac equation involving the Dirac-Coulomb (DC) Hamiltonian, 

Za 

H DC ip = E DC ip, H DC = ct ■ p + 13m . 

r 

Here, natural units (H = c = en = 1) are employed, and the familiar Dirac matrices are denoted 
by the symbols a and f3. The energy of an nS^M state (we use the usual spectroscopic notation 
for the quantum numbers) reads, when expanded up to sixth order in the parameter Za, 
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This is a power expansion in the parameter Za, where Z is the nuclear charge number and a is 
the fine-structure constant, and for Za < 1, it converges to the well-known exact Dirac-Coulomb 
eigenvalue [6]. 

On the other hand, let us suppose, hypothetically, that the electron were to carry no spin. 
Then, the equation would change to the bound-state equation for a Klein-Gordon particle, 

Z ex 

H KG ip = E KG i{), H KG = \fp 2 + m 2 . 

r 

In the expansion of an S-state energy levels in terms of Za, an irregularity develops for spinless 
particles, namely, a (Za) 5 term, and the (Za) e term carries a logarithm (see [7] for a detailed 
derivation) : 



(Za) 2 m (Za) A m(^ 3\ 8(Za) 5 m (Zafm 
E K g = m tt^ ^ — 1 ~~ tt~ M o q 1" 



7T 



2n 2 n 3 \ 8n J 3 ir n 3 n 3 

+ 0(Za) B . 



ln(Za) + ^C(3) 



2 . . n . . 1 5 5 29 

+ *(n) + 7e - ln(n) h 



2 ' n 3 n 2 16 n 3 12 



The expansion is nonanalytic (we denote by & the logarithmic derivative of the Gamma function, 
and 7e is Euler's constant). The occurrence of nonanalytic terms has been key not only to gene- 
ral bound-state calculations, but in particular also to Lamb shift calculations, which entail 
nonanalytic expansions in the electron-nucleus coupling strength Za in addition to power series 
in the quantum electrodynamic (QED) coupling a. A few anecdotes and curious stories are 
connected with the evaluation of higher-order logarithmic corrections to the Lamb shift [8, 9, 
10, 11]. The famous and well-known Bethe logarithm, by the way, is the nonlogarithmic (in Za) 
part of the energy shift in the order a(Za) A , and it is a subleading term following the leading- 
order effect which is of the functional form a(Za) 4 ln(Za). 

It does not take the additional complex structure of a Lamb shift calculation to necessitate the 
introduction of logarithms, as a simple model example based on an integral demonstrates [12], 







1-uj 2 ' [2 \pj 4 J I 16 \(3J 64 

Another typical functional form in the description of nature, characteristic of tunneling pheno- 
mena, is an exponential factor. Let us consider, following Oppenheimer [13], a hydrogen atom 
in an external electric field (with field strength \E\). The nonperturbative decay width due to 
tunneling is proportional to 



exp 



2(Za) 3 m 
3|e£| 



where \eE\ is the modulus of the electron's electric charge multiplied by the static electric field 
strength. 

We have by now encountered three functional forms which are typically necessary in order 
to describe expansions of physical quantities: these are simple powers, which are due to higher- 
order perturbations in some coupling parameter, logarithms due to some cutoff, and nonanalytic 
exponentials. The question may be asked as to whether phenomena exist whose description 
requires the use of all three mentioned functional forms within a single, generalized nonanalytic 
expansion? 
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The answer is affirmative, and indeed, for the description of energy levels of the double-well 
potential, it is known that we have to invoke a triple expansion in the quantities exp(— A/g) , ln(g) 
and powers of g in order to describe higher-order effects [14, 15] (here, g is a coupling parame- 
ter which is roughly proportional to the inverse distance of the two minima of the double-well 
potential). Other potentials, whose ground-state energy has a vanishing perturbative expan- 
sion to all orders (e.g., the Fokker-Planck potential), also can be described using generalized 
expansions [16]. The double-well and the Fokker-Planck Hamiltonians have stable, real eigenva- 
lues when acting on the Hilbert space of square- integrable wave functions (no complex resonance 
eigenvalues). An interesting class of recently studied potentials is PT-symmetric [17, 18, 19, 20]. 
Odd anharmonic oscillators for imaginary coupling fall into this class, but the double- well and 
the Fokker-Planck Hamiltonians do not. The purpose of this contribution is to assert that 
the concept of PT-symmetry is helpful as an auxiliary device in the study of odd anharmonic 
oscillators. 

In contrast to our recent investigation [21], we here focus on a few subtle issues associated 
with the formulation of the dispersion relation for odd anharmonic oscillators (Section 2), before 
giving a few new results for the cubic and quartic anharmonic oscillators in Section 3. In [21], 
by contrast, we focus on the sextic and septic oscillators. Conclusions are reserved for Section 4. 

2 Toward anharmonic oscillators 

Let us briefly recall why it is nontrivial to write dispersion relations for the energy levels of 
odd anharmonic oscillators. We consider as an example an odd perturbation of the form 7X 3 , 
with a coupling parameter 7, and we emphasize the differences to even anharmonic oscillators. 
Let us therefore investigate, as a function of the coupling parameter 7, the quartic and cubic 
potentials ^(7, x) = \x 2 + 7a; 4 and ^(7, x) = \x 2 + 7a; 3 . 

For the quartic potential and positive coupling 7 > 0, the spectrum of the Hamiltonian 
— \d 2 + n(7 > 0,a;), endowed with I? boundary conditions, consists of discrete energy levels 
which represent physically stable states with zero decay width. For 7 < 0, the potential ^(7, x) 
has a double-hump structure, and the particle can escape either to the left or to the right of 
the "middle valley" by tunneling (see Figs. 1(a) and 1(b)). So, when we change the sign of 
the coupling parameter, then "the physics of the potential changes drastically." We can then 
use the fact that, as a function of 7, the energy eigenvalues of the quartic oscillator have a 
branch cut along the negative real axis [1, 2, 3] and write a dispersion relation. It has been 
stressed in [22, 23] that the discontinuity of the energy levels is given exactly by the instanton 
configuration, and this fact has been widely used in the literature in the analysis of related 
problems in quantum physics and field theory. 

(Actually, when acting on L 2 , the negative-coupling quartic potential still possesses a real 
spectrum with discrete eigenvalues, but the analysis is highly nontrivial [24]. Indeed, the natural 
eigenenergies that are obtained from the real energies for positive coupling by analytic continu- 
ation as the complex argument of coupling parameter and of the boundary conditions, are just 
the complex resonance energies for which the dispersion relation holds.) 

Now let us investigate the odd potential 11(7, x) = x 2 /2 + 7X 3 . When 7 here changes sign, the 
physics of the potential does not change (see Figs. 1(c) and 1(d)): still, the particle can escape 
the "middle valley" by tunneling. Resonances occur. The question is whether we now have two 
branch cuts as a function of 7, one along the positive-7 axis and one for negative 7. Should 
we attempt to formulate a dispersion with integration along 7 £ (— 00, 0) and 7 £ (0,oo)? The 
answer is no. Rather, we should redefine the coupling in such a way that the PT-symmetry 
of the potential is used effectively. This means that the spectrum is real for purely imaginary 
coupling 7 = i (3 with real /3, and it is invariant under the transformation 7 = i [3 — > 7 = — i f3. 
In some sense, the case of the cubic potential for purely imaginary coupling is equivalent to the 
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(a) 



(b) 





(d) 



Figure 1. Plot of the potentials 14(7, x) 




• 4 and v(j, x) 




3 for positive and negative 



coupling 7 = ±1. 

quartic potential for positive coupling parameter, and the case of the cubic potential for positive 
coupling is equivalent to the quartic potential for negative coupling parameter. The key thus is 
to formulate the energy levels of the cubic as a function of g = 7 2 , not 7 itself [18]. 

3 Some results 

We here summarize a few results obtained recently [21] regarding the higher-order corrections for 
the energy levels of general even-order and odd-order anharmonic oscillators, using the quartic 
and cubic potentials as examples. Let us thus consider the two Hamiltonians, 



in the unstable region, i.e. for g < in the quartic case, and for g > in the cubic case. We 
assume both Hamiltonians to be endowed with boundary conditions for the resonance energies 
(which leads to a nonvanishing negative imaginary part for the resonance energy eigenvalues). 
Specifically, we denote the resonance eigenenergies by E n {g) for the quartic and e n (g) for the 
cubic, respectively. The quartic potential is plotted in the range g G (-2,-|) in Fig. 2(a), and 
the cubic potential is plotted in the range g 6 (§,2) in Fig. 2(b). 

Let us now investigate the instanton actions (see also Fig. 3). We write the classical Euclidean 
actions for the quartic and cubic, respectively, as 



H{g) 



Kg) 



2dx 2 



- + U(g,x), U{g,x) = -x 2 + gx 4 , 
+ V(g, x), V(g, x) = l -x 2 + ^gx 3 
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(a) (b) 

Figure 2. (a) The quartic potential U(g,x) = \x 2 + gx 4 is plotted in the parameter range g £ 
(—2.0,-0.5). As a function of the coupling g, the distance that the quantal particle has to tunnel 
before it escapes from the middle valley decreases as the coupling parameter increases. The decay width 
of the ground state is proportional to exp[l/(3g)] (for g < 0) and increases with the coupling. Higher- 
order corrections to this well-known result are indicated here (see equation (5)). (b) The cubic potential 
V(g,x) = \x 2 + i/gx 3 is plotted in the parameter range g £ (0.5,2.0). The decay width of the ground 
state is proportional to exp[— 2/(15g)] and increases with the coupling. 



and perform the following scale transformation x = x(t) = (—g) 1 ' 2 q(t) and y = y(t) = 
— g~ 1 ' 2 r(t) to arrive at 

S[q] = --Jtt (U 2 + W ~ 3 4 ) > 4r) =- g fdt (hr 2 + V ~ r 3 ) • 

Indeed, the width of the resonance is proportional to the exponential of minus the Euclidean 
action of the instanton configuration, which in turn is a solution to the classical equations of 
motion in the "inverted" potentials F(q) = q 4 — \q 2 and G(r) = r 3 — \r 2 . The instanton is 
given in Fig. 3. The instanton solutions read 

9(f) = 9ci(f) = ± [cosh(2t) + 1]" 1/2 (1) 



for the quartic and 



r{t) = r cl (t) = [cosh(i) + I]" 1 (2) 

for the cubic potential. Evaluating the instanton action, one obtains the leading-order results, 

1 2 

%cl] = — 5-1 s i r c\] = T^~- 
3g 15 g 

Observe that both instanton actions are positive in the relevant regions, where the potential is 
unstable (g < and g > 0, respectively). Consequently, the decay widths of the resonances of 
the quartic and cubic potential are proportional to exp[l/(3g>)] and exp[— 2/(15g)], respectively. 

In order to evaluate higher-order corrections are general formulas for oscillators of arbitrary 
degree, one needs dispersion relations. These reads for the quartic and the cubic, respectively, 

w ( \ j- 1 9 f° ^ Im£n(s + iO) ,„ , 

E n {g) =n+- / ds . — , (3a) 

2 ^J-oo s{s-g) 
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and [18] 

e n {g) = n+- + - ds r . (3b 

2 vr J s{s - g) 

One might ask if the integration for the cubic really stretches to s = +oo. The answer is 
affirmative: according to [25] , we may write the leading terms for the complex strong-coupling 
expansion for the first three resonances of the cubic as 

e ( 5 + iO) s =°° g 1 ^ 0.762851775 e^^ 5 , 
ei(fl + i0) 9 =°° 5 1/5 2.711079923e- i7r / 5 , 
e 2 (<? + i0) g =°° g 1 / 5 4.989240088 e- i7r / 5 . 

Here, we choose boundary conditions for the wave functions such as to generate resonance 
energies with a negative imaginary part, which are relevant for the dispersion integral (3b) as 
they are "attached" to values of the coupling with an infinitesimal positive imaginary part. 
Intuitively, we might assume that at least the second and the third resonance might disappear 
for very strong coupling g. This is because the classically forbidden region of the cubic potential 
which separates the fall-off region from the "middle valley" becomes smaller and smaller as the 
coupling increases, and indeed, the second excited level lies well above the relevant energy region 
in which tunneling would be necessary (see also Fig. 2(b)). However, this point of view does 
not hold: the resonance persists for arbitrarily large coupling, and the physical picture is that 
the "escape" of the probability density to infinity, which for the cubic happens in finite time, 
provides for a sufficient mechanism to induce a nonvanishing decay width of the resonance, 
even if the traditional tunneling picture is not applicable (similar considerations apply to the 
resonances in the Stark effect [26]). These considerations can be generalized to odd potentials of 
arbitrary order, and to arbitrary excited levels [21]. One result stemming from this generalization 
is given in the Appendix. 

Using generalized Bohr-Sommerfeld quantizations which are inspired by the treatment of 
double-well-like potentials [14, 15], one can formulate a general formalism [21] which allows to 
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write down higher-order formulas for the complex resonance energies. In contrast to [21], where 
we focused on the first few correction terms for the anharmonic oscillators of the third, sixth 
and seventh degree, here we would like to fully concentrate on the cubic and quartic oscillators 
and indicate the generalized nonanalytic expansion exclusively for the oscillators of the third 
and the forth degree. Specifically, we have for a resonance of the quartic, 



0) = J2 En > K - 



,K 



K=0 
oo 



2 2n+| 



j=i 

and for a general resonance of the cubic, 

oo 

tn{g > 0) = ^ e n,K9 K 

i 2 3 " 



vrn! (- g ) n +2 



exp 



1 

35 



J-i 



L=0 



~(4,n) K 



(4a) 



K=0 



K=0 

oo 

j=i 



r exp 



2 
1% 



J-i 

J L=0 



\^ -(3,n) if 

/ , ^j,L,i<y ' 



(4b) 



where the 3 are constant coefficients. (In contrast to [21], we here single out the perturbative 
contributions X}i?=o E n ,K9 K and Y^K=o e n,K9 K from the instanton effects, which are given by 
the terms with J = 1, . . . , oo.) 

Of particular phenomenological relevance is the term with J = 1 as it contains the pertur- 
bative corrections about the instanton configuration and is very important for comparison with 
numerically determined resonance eigenenergies of the systems. Without details, we only quote 
here [21] the results for the higher-order corrections to the ground state and to the first excited 
state of the quartic, which read 



lmE (g < 0) = - exp Q-^j 



11481557783 . 

9 + 



2 f 95 13259 . 
~ng\ + 24 9 ~ 1152 9 



8956043 , 

H 9 

82944 - 



[ 24 1152 82944 
4580883830443 5 12914334973382407 , 
iQiin9Q7fi 9 27518828544 9 



+ 



7962624 y ^ 191102976 zoxoozojw 

6938216714164463905 7 33483882026182043052421 , , ,, 

9 9 + 0(g 9 ) 

660451885056 126806761930752 



(5a) 



and 



32 



371 3371 2 
l + ^5-^5 2 + 



33467903 , 

q 

82944 



ImEUg < 0) = - exp ( - 

\3gJ V I 24" 1152" 82944 " 

73699079735 4 44874270156367 5 181465701024056263 6 

7962624 9 + 191102976 9 27518828544 9 ^ ' 

133606590325852428349 7 850916613482026035123397 8 9 1 

+ 660451885056 9 126806761930752 9 + ^ > ) ' 



660451885056 
and for the lowest two levels of the cubic, which are 



Ime (s > 0) 



exp - 



2 
15 S 



( 169 
~ ~16 9 ~ 512 



44507 2 
9 



128830328039451 5 1027625748709963623 , 
41943040 9 6710886400 " 



86071851 o 189244716209 4 

9 9 

y 2621440 y 



40960 



(6a) 
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and 



933142404651555165943 7 7583898146256325425743381 „ „ N 

9 9 + o(g ) 

107374182400 13743895347200 v ; 



8exp ( -jf^ j { 853 33349 2 395368511 8 



Imei(5 > 0) = . — <^ 1 - —5 + — g 



V^? 3 / 2 I 16 y 512 * 40960 y 

_ 1788829864593 4 _ 2121533029723423 5 

2621440 9 41943040 9 

27231734458812207783 6 37583589061337851179291 7 

6710886400 5 107374182400 5 
442771791224240926548268373 „ q ,l 
13743895347200 S + ° fe > j " (6b) 

Note that the higher-order terms for the ground state of the cubic, by virtue of the dispersion 
relation (3b), are in full agreement with the higher-order formulas given in [18]. Note also 
that both above results for the ground state could have been found by plain perturbation theory 
about the instanton configuration, but the results for the excited states are somewhat less trivial 
to obtain; they follow from the general formalism outlined in [21]. 



4 Conclusions 

Our generalized nonanalytic expansions (4a) and (4b) provide for an accurate description of 
resonance energies of the quartic and cubic anharmonic oscillators. These combine exponential 
factors, logarithms and power series in a systematic, but highly nonanalytic formula. Note that 
the term "resurgent functions" has been used in the mathematical literature [27, 28, 29] in 
order to describe such mathematical structures; we here attempt to denote them using a more 
descriptive, alternative name. 

In a general context, we conclude that a physical phenomenon sometimes cannot be described 
by a power series alone. We have to combine more than one functional form in order to write 
down a systematic, but not necessarily analytic expansion in order to describe the phenomenon in 
question. In the context of odd anharmonic oscillators, the generalized nonanalytic expansions 
which describe the energy levels in higher orders are intimately connected to the dispersion 
relations (3a) and (3b) which in turn profit from the PT-symmetry of the odd anharmonic 
oscillators for purely imaginary coupling. The 'PT-symmetry is used here as an indispensable, 
auxiliary device in our analysis (it is perhaps interesting to note that the use of 'PT-symmetry 
as an auxiliary device has recently helpful in a completely different context [30]). In our case, 
very large coefficients are obtained for, e.g., the perturbation about the instanton for the first 
excited state of the cubic (see equation (6b)). At a coupling of g = 0.01, the first correction term 
— 853g/16 halves the result for the decay width of the first excited state, and the higher-order 
terms are equally important. 

We have recently generalized the above treatment to higher-order corrections to anharmonic 
oscillators up to the tenth order. The oscillators of degree six and seven display very peculiar 
properties: for the sixth degree, some of the correction terms accidentally cancel, and for the 
septic potential, the corrections can be expressed in a natural way in terms of the golden ratio 
4> = (\/5 + l)/2. For the potential of the seventh degree, details are discussed in [21]. 

Let us conclude this article with two remarks regarding the necessity of using general non- 
analytic expansions to describe physical phenomena. First, the occurrence of the nonanalytic 
exponential terms is connected with the presence of branch cuts relevant to the description of 
physical quantities as a function of the coupling, as exemplified by the equations (3a) and (3b). 
Second, the presence of higher-order terms in the generalized expansions is due to our inability 
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to solve the eigenvalue equations exactly, or, in other words, to carry out WKB expansions in 
closed form to arbitrarily high order. These two facts, intertwined, give rise to the mathematical 
structures that we find here in equations (4a) and (4b). 



Appendix 



Using the dispersion relation (3b) and a generalization of the instanton configuration (2) to 
arbitrary odd oscillators, one may evaluate the decay width for a general state of an odd potential 
and general large-order ( "Bender- Wu" ) formulas for the large-order behavior of the perturbative 
coefficients of arbitrary excited levels for odd anharmonic oscillators. For a general perturbation 
of the form ^/gx M , with odd M > 3, with resonance energies e n M \g) ~ J2k e nK9 K '■> we 
obtain [21] in the limit K —* oo, 



(M) 

n.K 



(M - 2)r ((M - 2)K + n+\\ 

7r 3/2 n!2 2X+l-n 



B 



M 



M 



M-T M-2 



-(M-2)K-n- 



where B(x,y) = T(x)T(y)/T(x + y) is the Euler Beta function. 
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